Photon emission rate from atomic systems in the 

CSL model 

Stephen L. Adler 

Institute for Advanced Study 
Einstein Drive 
Princeton, NJ 08540 

Fethi M. Ramazanoglu 

Department of Physics 
Princeton University 
Princeton, NJ 08544 

Abstract 

We calculate the photon emission rate from a general atomic system in the 
mass-proportional CSL model. For an isolated charged particle emitting kilo- 
volt gamma rays, our results agree with those obtained by Fu. For a neutral 
atomic system, photon emission is strongly suppressed for photon wavelengths 
much larger than the atomic radius. However, for kilovolt gamma rays, Fu's 
result is modified by a structure factor that is of order unity, giving no rate sup- 
pression. Our calculation is readily generalized to the case of non-white noise, 
noise couplings that are not mass-proportional, and general (non-Gaussian) 
spatial correlation functions, and corresponding results are given. We briefly 
discuss the implications of our calculation for upper bounds on the CSL model 
parameters. 

1 Introduction 

Stochastic modifications of the Schrodinger equation, such as the continuous sponta- 
neous locahzation (CSL) model, solve the measurement problem in quantum theory 



2 



Adler and Ramazanoglu 



by giving an objective account of state vector reduction [1]. To assess the viability 
of these models, it is necessary to estimate lower and upper bounds on the stochastic 
model parameters, as surveyed in a recent paper of Adler [2]. An important upper 
bound on the stochastic rate parameter comes from a calculation by Fu [3] of the rate 
of noise-induced gamma radiation from free electrons, which he compares with the 
observed bound on 11 kilovolt gamma radiation from germanium. Adler suggested 
in [2] that in a neutral atomic system, radiation from protons, in the case of mass- 
proportional noise couplings, will largely cancel the radiation from electrons. Our aim 
in this paper is to check this assertion by a detailed calculation of stochastic noise- 
induced radiation in a general atomic system. We find that the asserted cancellation 
is present only for very long wavelength photons, whereas for the 11 kilovolt gamma 
rays figuring in Pu's bound, the radiation from protons somewhat enhances, rather 
than reducing, that from electrons. This result can be simply understood as the effect 
of inclusion of the space coordinate-dependent phase factor for the radiated wave. 

Thus for white noise, the upper bound on the CSL rate parameter is six orders of 
magnitude lower than estimated in [2], and hence is three orders of magnitude smaller 
that the lower bounds estimated in [2] from processes of latent image formation, 
assuming that latent image formation (and not subsequent development) corresponds 
to state vector reduction. Hence if the assumptions on which these lower bounds are 
based are correct, the white noise CSL model is disfavored. White noise is of course 
an idealization, and our calculation can be readily extended to the case of non-white 
noise. For non- white noise with a spectral cutoff below 11 kilovolts, there is no 11 
kilovolt gamma radiation, and so in this case the germanium experiment does not set 
a bound on the CSL model rate parameter, and there is no conflict with the lower 
bounds estimated in [2]. 

This paper is organized as follows. In Sec. 2 we outline the basic strategy of the 
calculation, which is to replace the real noise of the CSL model by an imaginary noise, 
that can be represented by a perturbation term in the Hamiltonian. We write down 
the general form of the Hamiltonian, and give the noise structure in the white-noise 
and non- white noise cases. 

In Sec. 3 we use standard atomic physics methods [4] to derive a master formula 
for the noise-induced photon radiation rate, in both the white noise and the non-white 
noise cases. In Sec. 4 we evaluate this formula for a single free electron, recovering the 
result of Fu [3] when his approximations are made. In Sec. 5 we evaluate the master 
formula for a hydrogenic atom, and in Sec. 6 for a general atomic system. In Sec. 
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7 we state the generalization of our results to a noise perturbation with general (not 
necessarily mass proportional) couplings to the particles, and with general spatial and 
time correlation functions. We conclude with a brief discussion of the implications of 
our calculation for CSL model phenomenology. 



2 General strategy, Hamiltonian, and noise struc- 



In the CSL model, the stochastic Schrodinger equation obeyed by the wave function 
ip takes the form dip = —{i/h)Hiljdt + N'lp + with H the usual Hamiltonian, 
with the noise term A/" real vahicd, and with the ellipsis ... representing additional 
nonlinear terms needed to preserve state vector normalization. A real valued choice 
for the noise term corresponds to an imaginary addition to the Hamiltonian, and is 
necessary to obtain a model that describes state vector reduction. An alternative 
stochastic Schrodinger equation can be written with an imaginary noise term, which 
docs not require additional nonlinear terms in the Schrodinger equation for norm 
preservation. This Schrodinger evolution docs not lead to state vector reduction, but 
for the case of white noise, it is a well known result that the noise average of the 
density matrix obeys the same evolution equation in the real and imaginary noise 
cases. Since the mean rate for noise induced transitions can be calculated from the 
noise averaged density matrix, this implies that one can use the imaginary noise 
equation to calculate the mean rate for such transitions. Hence, to leading order, one 
can represent the noise perturbation as a self-adjoint perturbation on the Hamiltonian 
H, and use standard second order perturbation theory to evaluate its effects. 

The usual justification for the use of imaginary noise is based on a calculation 
of the density matrix evolution in the real and imaginary noise cases using the Ito 
calculus, which as already noted, assumes white noise. Adler and Bassi [5] have 
recently shown, however, that in the case of non-white Gaussian noise, the noise 
averaged density matrix evolutions are still the same for the real and imaginary noise 
cases, through second order in the noise term. Hence, in the second order perturbation 
calculations of this paper, we can use an imaginary noise term to calculate the effects 
of non-white noise as well as white noise. 

We will thus be considering a Hamiltonian of the form 



ture 



H = Ho + Hem + 



n 
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with Hq the atomic system Hamiltonian, Hem the electromagnetic perturbation de- 
scribing photon emission, and iJ„ the perturbation describing the noise. For a system 
of particles of charges ej and masses rrij, the electromagnetic perturbation is 



Hem = y2^A{^j)-^x,+0{A') , (2) 

ITijC 
j=l J 

with the electromagnetic potential, for field quantization in a cubical box of size L, 
given by 

p V ^ 

where Up = pc, and where the numerical value of a unit unrationalized charge e is 
e^/{hc) ~ 1/137.04. Since we are only interested in the matrix element for emitting 
a single photon of wave number p, we pull this term out from Eq. [3] and, separating 
off the time dependence, write the electromagnetic perturbation as 



^ y ojpL-^ ^ nijC 

where V-,- is an abbreviation for V^;^.. 

In the CSL model with mass-proportional couplings, the noise perturbation can 
be written as 

V[z,{x})= 2^mjg{z-Xj) . 

j 

Here g{x) is a spatial correlation function, conventionally taken as the Gaussian 



with = Tc, and with conventionally taken as 10~^ cm. In the case of white 
noise, dWt is an Ito calculus differential that obeys 



dWt{x)dWt{y) = -fdt6\x-y) 



(7) 
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with 7 the noise strength parameter. The corresponding formula for the case of 
non-white noise is 

POO 

/ du^{uj)e-'^^'-''^S%x-y} , (8) 

' — oo 

with E[...] denoting the expectation or average over the noise. When 7(0;) is a constant 
7, Eq. [H] reduces, on integration over t', to Eq. [3 



E 



dWAx)dWH( 



dt 



dt' 



3 Master equation for the radiation rate 

According to Eqs. H] and [5], the total perturbation on the atomic Hamiltonian Hq is 



V{t) = j dh^^^V{z, {x}) + e'"''*>VP({f }) . (9) 

Expanding the transition amplitude in a perturbation series following the methods of 
[4], we get 

{f\Uj{t,0)\z) = l + T}p + TlP + ... , 



= -1 fds f duY,{f\Vi{s)\k) {k\Vj{u) |z) 

!>' Jo Jo u 



• [' ds f in r dE ei('=,--).e*<--.)« Y. 

Jo Jo J-00 u 

(10) 



2TTh Jo Jo J-00 ^Ei + ir]-Ek 



(2) 

where in the first line of the formula for Tj- , Vj denotes the interaction picture pertur- 
bation, and in the second line Vfk and Vki denote matrix elements of the Schrodinger 
picture perturbation. To calculate the noise induced radiation, we are only interested 
in the terms in Eq.[TO]that are bilinear in the electromagnetic and noise perturbations, 
so on substituting Eq. [9] and dropping irrelevant terms, we get 

poo 



00 



Jo J ds Jq ^ E + IT] — Ek \ > 



^ E + ir)-Ek 
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Taking the squared modulus of Eq. [TTl averaging over the noise, and using the for- 
mulas for representations of the Dirac delta function given in [4], in the large time 
limit we obtain in the white noise case, 

with the corresponding equation in the non-white noise case taking the form 

(13) 

Equations [12] and [13] are the master equations from which we shall calculate the noise 
induced radiation rate, by substituting the matrix elements of V and appropriate 
to the various cases of interest. 



+ 



Ei - huop + 17] - Ek Ef + huop + ir] - Ek 



(12) 



^ Ei - hup + ir] 



Ek Ef + hUp + ir] — Ek 



4 Free electron: repeating Fu's calculation 

As a first application of Eq. [12], and a check, let us repeat the calculation of Fu [3] 
for the case of a single free electron. Assuming that the electron is initially at rest, 
the initial, final, and intermediate state electron wave functions are 

i^i = -7= , ^/ = -7= , i'k = -7= ■ (14) 



From Eqs. 4 and 5, as specialized to a single particle of charge e (with e^/ {he) ^ 1/137) 
and mass m, the needed matrix elements are 



27r^c eh^ (15) 



and 
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Substituting these into Eq. [121 we get for the noise averaged squared matrix 
element 



mNL^\l pL^ c ^ -^(^p2 _^2p- q) + hcp + ir] 

Fu notes that when the photon momentum p is much larger than the inverse 
correlation length l/vc, the Gaussian factor in Eq. [T7] forces the electron and photon 
to emerge nearly back to back, that is, q ~ —p. As a result 

^ ^ (18) 
-£^{p' + 2p.q) ' 

which is of order 100 for hpc =11 keV. Thus one can to a good approximation keep 

only the term hep in the denominator of Eq. [T71 which then simplifies to 

Integrating over phase space for the electron and photon, summing over photon po- 
larizations, and dividing by the elapsed time, we get for the radiated power per unit 
photon momentum space volume and per unit time. 



• (17) 



(20) 



Carrying out the integrals and polarization sum, and replacing the noise parameter 7 
by a new parameter A defined by 7 = Svr^/^rj^A, we get finally for the power radiation 
rate 

dp irr^mj^p 

This is in agreement with the result obtained by Fu [3], when our unrationalized 
charge squared is replaced by e^/(47r), corresponding to Fu's use of a rationalized 
charge convention. 



5 Hydrogenic atom 

We consider next a hydrogenic atom, with oppositely charged particles of masses mi 
and 7712. Equation 5 for V{z, {x}) now takes the form 

V{z,{x})=-—M{z,{x}) , (22) 

M{z,{x}) =mig{z- xi) + m2g{z- X2) ■ (23) 
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Introducing the center of mass coordinate X, total mass M, relative coordinate x, 
and reduced mass n, by 



- mi ^ m2 ^ 

M— 1711+1712, H 



X ^ Xi - X2 , 
17111712 



M 



(24) 
(25) 



we can use the fact that the Bohr radius Oq is much smaller than Vc to approximate 
M.{z, {x}) as follows, 



M{z,{x}) ^mig{z-xi) +m2g{z-X2) 

17111712 ^ .2 



Mg{z-X) + 
Mg{z-X) , 



2M 



ix-V,yg{z-X) 



(26) 



givmg 



27:hc ieTi 



-ip-xi 



mi 



Vi e 

m2 



-tp-X2 



V{z, {x}) 



hM 

rriN 



g{z-X) . 



The initial, final, and itermediate state atomic wave functions are now 



1 



Jk-X 



where we use carets to denote the labels of hydrogenic internal states. Defining 



(27) 
(28) 

(29) 



0{k) = -[e- 



■ mi -* -* 



ep-k+ (—e-''^P-^ + — e^^^^- ) ■ V, , (30) 
\mi i7i2 ' 



we find that the matrix elements entering the master formula are 

O{0) 



2-11 he ieh I , 
— < k 



2TThc ieh ^ 
pL^ e 



0{k) 



k) d 



k—p—q 1 



and 



hM 



hM 



ik-z-\k'^rl 



(31) 
(32) 

(33) 
(34) 
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Then without any further approximation we find 

^ 27rhc e^h^ 



-ft f hM Y 



P 



f o{6) 



+ 



/ 



0{p + ^ 



2M 



hep 



2M ^ "/ 



(35) 



with Efi = — El the internal energy difference between the final and initial atomic 
states. The radiated power, per unit photon momentum space volume and per unit 
time, now requires a sum over the final internal atomic state /, and is given by 

' ' " ; " [ (36) 



sian e 



Note that when p + q = 0, the two terms in Eq. 37 cancel. Since the Gaus- 
constrains |p + ^ to be not much larger than l/vc, we can make this 



cancellation explicit by expanding in the small parameter 

fi^p-(p + q) _%^p-{p + q) 

which keeping the leading two terms, and writing p ■ x = pz, gives 



(37) 



dP 



dp c 



MVe^X 



rriN 



f 



f 



,2fe4 



+ 



p^h 



f 



1 



-e M 



pz 



d 



mi 777.2 / dx 

For small p, this expression can be further simplified to 



(38) 



dP_ 
dp 



p 



777JV / Trr^ 



1 



f 



+ 



M^E%t,^ 



f 



d_ 
dx 



.h^ 1 e^Av- 

2P T — 

c m,if 7rr^ 



(39) 



Eh 



where we have used the dipole approximation formula 



/ 



dx 



f 



x 



(40) 
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which shows that the two terms in Eq. [39] make equal contributions. The sum in 
Eq. [39] has been evaluated in closed form by Dalgarno and Kingston [6], with the 
result 



/ 



/ 



2^6 



^2 



43 fi^d, 



(41) 



giving an explicit expression for the small p radiation rate. 

However, for 11 kilovolt photons, the small p approximation does not apply, and 
instead we can simplify the formulas by making the approximation Dq ^ hep, as was 
done by Fu in his calculation. The radiation rate then becomes 

2 9^ 



dP 
dp 



I h f M 



p c \mj^ 



7rr2 ^ 1 M2 V 



f 



mo 



e M 



pz 



+ 



f 



nil 



M -| e' M 

m2 



1, 



pz 



d_ 

dx 



(42) 



which, using the completeness of the hydrogen spectrum, can be simplified to 

2 9^ 



dp 

dp 



I h ( M 



p C \mj^ 



+ 



f 



M2 



2 — 2 cospz 



d 



dx \mi 



2 \ d 
cospz — 

17111712 J OX 



(43) 



The ratio of the second term to the first can be shown to be of order (e^/Zic)^, so 
the second term can be neglected. Evaluating the first term using the hydrogen atom 
ground state wave function, we find the final result for high p to be 



dP 
dp 



1 - 



1 



1 h 1 e^A 



p m\ T^rl 



(44) 



For small p this expression is suppressed with respect to the rate calculated by Fu, 
but for large p it approaches twice Fu's rate, because when the photon wave length 
is much smaller than the atomic radius, the electron and proton radiation rates add 
incoherently. For 11 kilovolt gamma radiation from hydrogen, the rate given by Eq. [H] 
is about 1.8 times the rate for a free electron. The structure of the first term in Eqs. [12] 
and [13] can be readily understood in terms of the phase factor that appears in the 
formula for the radiation rate of a distributed charge system, as in Eqs. (13-33) and 
(13-37) of the text of Panofsky and Phillips [7]. 
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6 Many-body system 

We turn next to a general n particle atomic system, for which the electromagnetic 
and noise perturbations are given by Eqs. 4 and 5, with the sum over j extending 
from 1 to n. In order to take account of overall momentum conservation, we separate 
the coordinates of the particles into a center of mass coordinate X and internal 
coordinates ,i = 1, ...,n — 1, hj writing 



n 

Xi = + X i = I, n - 1 

n—l 

I 

Xn = X 



n—l 



m„ ^ (45) 



-^i ^ ] 



n _, 

rrijXj 



M 



M 

In the following equations, Vj denotes the partial derivative with respect to the 
original coordinate Xj, not the derivative with respect to the internal coordinate 
Straightforward calculations show that the commutator of this partial derivative with 
an internal coordinate is given by 

[a-V„6-0] = a-6('^.,-^) , (46) 
and also that the Jacobian J of the transformation of Eq. HS] is given by 

J = ^gi:4^ = (-1)"- f 1 + . (47) 



9(xei...en-i 



Moreover, the kinetic term of the unperturbed hamiltonian is separated by the trans- 
formation of Eq. Us] into a center of mass part and an internal part, 

y^ =— yvj , (48) 



2mi 2M ^ 2mi 2M , 

i=l i=l \j=l 



12 



Adler and Ramazanoglu 



so that we know that wave functions are of the factorized form 



1 



Ak-X 



(49) 



Using the center of mass transformation and the factorized wave functions, the 
noise and radiation matrix elements needed for the master formula of Eq. 12 are 
calculated to be 



pL^ c \ I ^ rrii 



S 



k+p 



,ep ■ k 



m 



and 



Vki{z) 



Vfk{?) 



J{k-q)-z-Uk~q}'^rl 



/IE' 



-i{k-q)-£,j 



(50) 
(51) 

(52) 
(53) 



We now simplify Eq. 12 by making the approximation that the photon energy 
hujp is much larger than both the internal energy differences and the center of mass 
recoil energy, that is, that hup is much larger than Ei — Ek and Ef — Ek. With this 
approximation (which is analogous to the approximation made by Fu and also made 
in Eqs. 42-44 of our hydrogen atom calculation), Eq. 12 simplifies to 

2 



E 



(54) 

k ■ ~-^P 

Substituting Eqs. 50-53 into Eq. EH summing over the final state / by the analog 
of Eq. 36, and using completeness twice together with algebraic simplification using 
Eq. Hni we get for the power radiated 



dP _ 2-i h 1 r dVtp 
dp (2vr)^ m^c^ p J Att 



d w e 



{w ■ py 



(55) 



Note that the internal integration to be used in evaluating the matrix element in this 
formula includes the Jacobian J of Eq. HTl and so is 

' \, . (56) 



n — 1 » 

ji n / d% 
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To check that Eq. [55] reproduces the result of the first term of Eq. 43 for the 
hydrogen atom, we note first that for a two particle system one has xi = X + ,X2 = 
X + ^2, and so X = Xi — X2 = ^1 — ^2, which by Eq. reduces to x = ^i(l + mi/m2). 
Hence |J|(i^^i = (1 + mi/m2)^d^C,i = d^x, so the internal integration involves the 
conventional internal coordinate used for the hydrogen atom. The expansion in the 
small parameter of Eq. 37 is equivalent, in the many-body context, to setting w = 
in A/" in Eq. [55l an approximation that permits the integration over w to be easily 
done, yielding our previous formula for the hydrogen atom radiated power. 

One can also apply Eq. to the case of a crystal lattice. Again making the 
approximation of neglecting w in Af, that is, taking Vc to be large, we define 



(57) 



cell 



We then find that the matrix element appearing in Eq. takes the form (with (...) 



denoting an expectation in the initial state 



and with Li a lattice displacement). 



(iA/-p> = A'„„((i/r)-i{/>r) + 



-ip-Li 



l(/)P 



(5^ 



= iv.e.(i/-(/)r) , 

since the second term on the first line of Eq. [58] grows more slowly than Nceii for 
generic values of p. Hence as long as the variance of / over a unit cell is nonzero, the 
radiated power scales as the size of the crystal lattice (at least for lattice dimensions 
smaller than Vr). 



7 Generalizations and discussion 

Several generalizations of the formulas given above can be easily derived. First of all, 
if the noise Hamliltonian of Eq. 5 involves general couplings Qi that may differ from 
the masses mj, so that 



h (59) 

V{z,{x})= y29j9{z~Xj) , 

3 

then in Af in Eq. [55] one replaces Cj by ejgj/mj. Secondly, our calculation, in the non- 
white noise case, can be viewed as calculating the radiation produced by a random 
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gravitational potential 

Vgrav{x,t) ='^mi(f){Xi,t) , (60) 

i 

with {4>)av = and with the correlation function 

X, t)(f){x', t'))Av = / du;^{uj)e-'-^''''^G{x - x') , (61) 

\mNj 271 

G{x-x') = I d^zg{x- z)g{z-x') . (62) 



Since for the Gaussian g of Eq. 6 one has 

d'^xe^^-^gix) =e-^^'^' , (63) 

/.».e.^^-G(f) ^l^.,'-^l^y,i,-my> - e-"'' . (64) 

for a general G{x) in Eq. [61] one simply replaces e~"'^^= in the radiated power expres- 
sions by 

G[w] = [ c/=^xe™-"^G'(f) . (65) 



Finally, for a more general non-white noise that does not have a time-translation 
invariant correlation function, so that Eq. 8 is replaced by 



E 



dWAx) dWA 



A{t,t')6%x-y) , (66) 



dt dt' 

the master formula in the non-white noise case is modified by replacing 

t^{u, + El^) (67) 

by 

ds [ rfMA(s,M)e*(^-")["''+^^l . (68) 
Jo 

The most general case, in which the correlation function of Eq. 61 does not factorize 
into a temporal correlation times a spatial correlation, can be obtained by combining 
results from Eqs. 61-68. 

To conclude, we consider the implications of our results for CSL model phe- 
nomenology. Since we have seen that for a hydrogenic or a general atomic system 
emitting kilovolt gamma rays, charge neutrality does not imply a corresponding can- 
cellation in the radiation rate, the estimates of Fu [3] must be taken as giving the 
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best upper bounds on the CSL parameter A (defined following Eq. 20) for the white 
noise case. Including [2] a factor of An correction to Fu's evaluation of the electric 
charge squared e^, as well as [8] a factor of roughly 4 increase in the experimental 
rate limit subsequent to the value used by Fu, Fu's calculation implies the bound 
A < 7 X 10~^^s~^, which is ~ 3 x 10^ larger than the standard CSL model value of 
A = 2.2 X 10~^^s~^. As we noted in Sec. 1, this upper bound is several orders of 
magnitude below the lower bound on A set by postulating that latent image forma- 
tion (as opposed to image development) should correspond to state vector reduction. 
Although increasing Tc to 10~^ cm decreases the 11 kilovolt photon radiation rate, 
and so increases the corresponding upper bound on A, by two orders of magnitude, 
as discussed in [2] this increase in Tc also increases the latent image formation lower 
bound on A by one to two orders of magnitude, and so does not eliminate the potential 
discrepancy. 

By contrast, in the non-white noise case there is not necessarily a conflict, since 
the relevant radiation rate involves the noise spectral coefficient '~f{uj) at a frequency 
of at least that of the emitted gamma ray, of order lO^^s"^ In fact, in their review 
[1], Bassi and Ghirardi suggest a cutoff in the noise frequency spectrum of order 
c/tc ~ 10^^ s~^, which would be more than sufficient. Even a much lower frequency 
cutoff would suffice to explain reduction in typical measurements with measurement 
times of order a nanosecond or longer; for example, a cutoff of order lO^^s^^ would 
be more than adequate. This would correspond to an energy cutoff of order 10^'' eV, 
or a noise temperature of order 1 degree K. So possibly even a non-white cosmic relic 
background noise field, with suitable correlator structure, coupling as a real-valued 
noise term A/" in the Schrodinger equation for dip, could explain state vector reduction 
in measurement situations, without coming close to violating the upper bound set by 
Fu's calculation. 
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9 Added Note 

The use of the term "power" and the symbol P in Eqs. (20), (21), (36), (38), (39), 
(42), (43), and (55) was inadvertent; we should have used the term "rate" and the 
conventional symbol F. These formulas all give the photon radiation rate, and do not 
include the energy per photon factor hep needed to convert them to radiated power. 
We wish to thank Angelo Bassi for pointing this out to us. 
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